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Abstract 

In this note we briefly review the recent studies of dark matter in the MSSM and its singlet 
extensions: the NMSSM, the nMSSM, and the general singlet extension. Under the new detection 
results of CDMS II, XENON, CoGeNT and PAMELA, we find that (i) the latest detection results 
can exclude a large part of the parameter space which allowed by current collider constraints in 
these models. The future SuperCDMS and XENON can cover most of the allowed parameter 
space; (ii) the singlet sector will decouple from the MSSM-like sector in the NMSSM, however, 
singlet sector makes the nMSSM quite different from the MSSM; (iii) the NMSSM can allow light 
dark matter at several GeV exists. Light CP-even or CP-odd Higgs boson must be present so as 
to satisfy the measured dark matter relic density. In case of the presence of a light CP-even Higgs 
boson, the light neutralino dark matter can explain the CoGeNT and DAMA/LIBRA results; (iv) 
the general singlet extension of the MSSM gives a perfect explanation for both the relic density 
and the PAMELA result through the Sommerfeld-enhanced annihilation. Higgs decays in different 
scenario are also studied. 

PACS numbers: 14.80.Ly,11.30.Pb,95.35.+d 
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I. INTRODUCTION 



Although there are many theoretical or aesthetical arguments for the necessity of TeV- 
scale new physics, the most convincing evidence is from the WMAP (Wilkinson Microwave 
Anisotropy Probe) observation of the cosmic cold dark matter, which naturally indicates 
the existence of WIMPs (Weakly Interacting Massive Particle) beyond the prediction of the 
Standard Model (SM). By contrast, the neutrino oscillations may rather imply trivial new 
physics (plainly adding right-handed neutrinos to the SM) or new physics at some very high 
see-saw scale unaccessible to any foreseeable colliders. Therefore, the TeV-scale new physics 
to be unraveled at the Large Hadron Collider (LHC) is the most likely related to the WIMP 
dark matter. 

If WIMP dark matter is chosen by nature, it will give a strong argument for low-energy 
supersymmetry (SUSY) with R-parity which can give a good candidate. Nevertheless, SUSY 
is motivated for solving the hierarchy problem elegantly. It can also solve other puzzles of 
the SM, such as the 3<r deviation of the muon anomalous magnetic moment from the SM 
prediction. In the framework of SUSY, the most intensively studied model is the minimal 
supersymmetric standard model (MSSM) [1], which is the most economical realization of 
SUSY. However, this model suffers from the /i-problem. The /z-parameter is the only di- 
mensional parameter in the SUSY conserving sector. From a top down view, one would 
expect the fi to be either zero or at the Planck scale. But in the MSSM, the relation of the 
electro-weak (EW) scale soft parameters (m^, ml) [2] 



1 M 2 = m 2 d -mlt&n z 
2 z tan 2 /3-l 



M 2 Z = \ ' -fi 2 , (1) 



makes that fi must be at the EW scale, while LEP constraints on the chargino mass require 
fi to be non-zero [3] . A simple solution is to promote fi to a dynamical field in extensions of 
the MSSM that contain an additional singlet superfield S which does not interact with the 
MSSM fields other than the two Higgs doublets. An effective /i can be reasonably got at EW 
scale when S denotes the vacuum expectation value (VEV) of the singlet field. Among these 
extension models the next-to-minimal supersymmetric model (NMSSM) [4] and the nearly 
minimal supersymmetric model (nMSSM) [5, 6] caused much attention recently. Note that 
the little hierarchy problem which is also a trouble of the MSSM is relaxed greatly in the 
NMSSM. 
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If introduce a singlet superfield to the MSSM, the Higgs sector will have one more CP 
even component and one more CP odd component, and the neutralino sector will have one 
more singlino component. These singlet multiplets compose a "singlet sector" of the MSSM. 
It can make the phenomenologies of SUSY dark matter and Higgs different from the MSSM. 
More and more precision results of dark matter detection give us an opportunity to test if this 
singlet sector really exists. For example, experiments for the underground direct detection 
of cold dark matter x have recently made significant progress. While the null observation of 
X in the CDMS and XENON100 experiments has set rather tight upper limits on the spin- 
independent (SI) cross section of x-nucleon scattering [7, 8]. The CoGeNT experiment [9] 
reported an excess which cannot be explained by any known background sources but seems 
to be consistent with the signal of a light x with mass around 10 GeV and scattering rate 
(1-2) x 1CT 40 cm 2 . Intriguingly, this range of mass and scattering rate are compatible with 
dark matter explanation for both the DAMA/LIBRA data and the preliminary CRESST 
data [10]. Though CoGeNT result is not consistent with the CDMS or XENON results, it 
implies that the mass of dark matter can range a very long scope at EW scale, that is from 
a few GeV to several TeV. The indirect detection PAMELA also observed an excess of the 
cosmic ray positron in the energy range 10-100 GeV [11] which may be explained by dark 
matter. 

In this paper, We will give a short review on the difference between the MSSM and the 
MSSM with a singlet sector under the constraints of new dark matter detection results. 
As the Higgs hunting on colliders has delicate relation with dark matter detections, the 
implication on Higgs searching is also reviewed. The content is based on our previous work 
[12-14]. the paper is organized as following, in sec. II, we will give a short review on the 
structures of the MSSM, the NMSSM and the nMSSM. In sec. Ill we will give a comparison 
on the models under the constraints of CDMS, XENON, and CoGeNT. In sec. IV, a general 
singlet extension of the MSSM is discussed, and a summary is given in sec. V. 

II. THE MSSM AND ITS SINGLET EXTENSIONS 

As the economical realization of supersymmetry, the MSSM has the minimal content of 
particles, while the NMSSM and the nMSSM extend the MSSM by only adding one singlet 
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Higgs superfield S. The difference between these models is reflected in their superpotential: 

W F + fiH u -H d , (2) 

(3) 
(4) 



Wmssm 
Hnmssm 

WnMSSM 



W F + \H u -H d S + -KS\ 



W F + \H U ■ H d S + Z F M*S, 



where W F = Y U Q ■ H U U - Y d Q ■ H d D - Y e L ■ H d E with Q, U and D being the squark 
superfields, and L and E being the slepton superfields. H u and H d are the Higgs doublet 
superfields, A, k and £p are dimensionless coefficients, and /x and M n are parameters with 
mass dimension. Note that there is no explicit /i-term in the NMSSM or the nMSSM, and 
an effective /i-parameter (denoted as /i e ff) can be generated when the scalar component (S) 
of S develops a VEV. Also note that the nMSSM differs from the NMSSM in the last term 
with the trilinear singlet term kS 3 of the NMSSM replaced by the tadpole term ^ F M^S. 
As pointed out in Ref. [5], such a tadpole term can be generated at a high loop level and 
naturally be of the SUSY breaking scale. The advantage of such replacement is the nMSSM 
has no discrete symmetry thus free of the domain wall problem which the NMSSM suffers 
from. 

Corresponding to the superpotential, the Higgs soft terms in the scalar potentials are also 
different between the three models (the soft terms for gauginos and sfermions are the same 
thus not listed here) 
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After the scalar fields H u ,H d and S develop their VEVs v u , v d and s respectively, they can 
be expanded as 

//.. 
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^75 (v d + (p d + i<p d ) 
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(v u + <Pu + i<Pu) 



— (s + a + it). (8) 



The mass eigenstates can be obtained by unitary rotations 
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where ^1,2,3 and are respectively the CP-even and CP-odd neutral Higgs bosons, G° 
and G + are Goldstone bosons, and H + is the charged Higgs boson. Including the scalar 
part of the singlet sector in the NMSSM and the nMSSM leads to a pair of charged Higgs 
bosons, three CP-even and two CP-odd neutral Higgs bosons. In the MSSM, we only have 
two CP-even and one CP-odd neutral Higgs bosons in addition to a pair of charged Higgs 
bosons. 

The MSSM predicts four neutralinos xi (i — 1 ; 2, 3, 4), i.e. the mixture of neutral gauginos 
(bino A' and neutral wino A 3 ) and neutral higgsinos (V'h,,) V'hJj wnne the NMSSM and the 
nMSSM predict one more neutralino corresponding to the singlino ips from the fermion part 
of singlet sector. In the basis (— iA^ — iA 3 2, ^s) (for the MSSM ips is absent) the 

neutralino mass matrix is given by 

/ Mi m z swSb —mzSwCb n ^ 

M 2 -m z c w s b m z c w c b 

m z swSb —m z swSb — /i 

-fJL 



-m z s w c b -m z c w c b 




— Xvcb —Xvsb 



\ 



2^ 







—Xvcb 
—Xvsb 

for the NMSSM 
for the nMSSM 



(10) 



/ 



where Mi and M 2 are respectively U(l) and SU(2) soft gaugino mass parameters, sy/ = 
sin9 w , cw = cos9 w , s b = sin/? and Cb = cos/3 with tan/3 = v u jv^. The lightest neutralino 
Xi is assumed to be the lightest supersymmetric particle (LSP), serving as the SUSY dark 
matter particle. It is composed by 



X\ = N u (-iX') + N u (-tX 3 ) + N^° Hu + N u ^ Hd + N^ s , 



(11) 



where is the unitary matrix (N15 is zero for the MSSM) to diagonalize the mass matrix 
in Eq. (10). For the mass matrices above we should note that the following two points 

1. For a moderate value of k, the neutralino sector of the NMSSM can go back to the 
MSSM when A approaches to zero. This is because in such case the singlino component 
will become super heavy and decouple from EW scale. The singlet scalar will not mix 
with the two Higgs doublet, then the NMSSM will be almost the same as the MSSM 
at EW scale. 
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2. Since the ips^s element of Eq. (10) is zero in the nMSSM, the singlino will not 
decouple when A approaches to zero. In fact, in the nMSSM the mass of the LSP can 
be written as 

_ 2fi cS \ 2 (v 2 u + v 2 ) tan/3 
m *°~ 2^ + \ 2 (v 2 u + v 2 )t^ 2 p + Y 1 } 

This formula shows that to get a heavy xl, we need a large A, a small tan/3 as well as 
a moderate fj, e g. 

The chargino sector of these three models is the same except that in the NMSSM/nMSSM 
the parameter \x is replaced by ii e s. The charginos xt,2 ( m xf — m x0 are ^ e m i x t ure of 
charged Higgsinos d and winos = (A 1 ± A 2 )/^, whose mass matrix in the basis of 
{-i\ ± ,^ u d ) is given by 

(M 2 \/2m w s b 

So the chargino xt can De wino-dominant (when M 2 is much smaller than /i) or higgsino- 
dominant (when \i is much smaller than M 2 ). Since the composing property (wino-like, 
bino-like, higgsino-like or singlino-like) of the LSP and the chargino xt is ver Y important in 
SUSY phenomenologies, we will show such a property in our following study. 




III. COMPARISON WITH THE MSSM AND THE MSSM WITH A SINGLET 
SECTOR 

A. In light of CDMS II and XENON 

First let's see the MSSM, the NMSSM and the nMSSM under the constraints of results 
of CDMS II and XENON100. As both current and future limits of x-nucleon of CDMS and 
XENON are similar to each other, we will show only one of them. Nevertheless, as a good 
substitute of the SM, SUSY model must satisfy all the results of current collider and detector 
measurements. In our study we consider the following experimental constraints: [15] (1) we 
require Xi t° account for dark matter relic density 0.105 < Vth 2 < 0.119; (2) we require 
the SUSY contribution to explain the deviation of the muon a M , i.e., a^ xp — a^ M = (25.5 ± 
8.0) x 10~ 10 , at la level; (3) the LEP-I bound on the invisible Z-decay, T(Z ->• X iXi) < 1-76 

6 



MeV, and the LEP-II upper bound on a(e + e~ — > which is 5 x 10~ 2 pb for i > 1, as 

well as the lower mass bounds on sparticles from direct searches at LEP and the Tevatron; 
(4) the constraints from the direct search for Higgs bosons at LEP-II, including the decay 
modes h — > h\h\,a\a\ — > 4/, which limit all possible channels for the production of the 
Higgs bosons; (5) the constraints from B physics observable such as B — > X s ^, B s — > 
B + — > t + u, T 7CI1, the ai^rjb mixing and the mass difference AM^ and AM S ; (6) the 
constraints from the precision EW observable such as pi cpt , sin 2 ^^*, m w and R^, (7) the 
constraints from the decay T — > 7/^, and the Tevatron search for a light Higgs boson via 4/i 
and 2/x2r signals [16]. The constraints (1-5) have been encoded in the package NMSSMTools 
[17]. We use this package in our calculation and extend it by adding the constraints (6, 7). 
As pointed out in Ref. [16], the constraints (7) are important for a light Higgs boson. In 
addition to the above experimental limits, we also consider the constraint from the stability 
of the Higgs potential, which requires that the physical vacuum of the Higgs potential with 
non- vanishing VEVs of Higgs scalars should be lower than any local minima. 

For the calculation of cross section of x-nucleon scattering, we use the formulas in Ref. 
[18, 19] for the MSSM and extend them to the NMSSM/nMSSM. It is sufficient to consider 
only the SI interactions between Xi and nucleon (denoted by f p for proton and f n for neutron 
[19]) in the calculation. The leading order of these interactions are induced by exchanging 
the SM-like Higgs boson at tree level. For moderately light Higgs bosons, f p is approximated 
by [19] (similarly for/ n ) 

fp- Y] — m pfT, + 7^fT G J2 ~ m P ' ( 14 ) 

q=u,a,s H q=c,b,t H 

where denotes the fraction of m p (proton mass) from the light quark q while fx G = 
1 — ^ u d s is the heavy quark contribution through gluon exchange, is the coefficient 
of the effective scalar operator. The x°-nucleus scattering rate is then given by [19] 




where mx is the mass of target nucleus and n p {n n ) is the number of proton (neutron) 
in the target nucleus. In our numerical calculations we take = 0.023, = 0.034, 
= 0.019, = 0.041 and = = 0.38. Note that the scattering rate is very 
sensitive to the value of f Ts [21]. Recent lattice simulation [22] gave a much smaller value 
of fr s (0.020), it reduces the scattering rate significantly which can be seen in Ref. [23]. 
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Considering all the constraints listed above, we scan over the parameters in the following 
ranges 

100 GeV < (Mg, M h m A , /i) < 1 TeV, 
50 GeV < Mi < 1 TeV, 1 < tan/3 < 40, 

(|A|,|k|)<0.7, |A K |<lTeV, (16) 

where Mg and M^ are the universal soft mass parameters of the first two generations of 
squarks and the three generations of sleptons respectively. To reduce the number of the 
relevant soft parameters, we worked in the so-called m™ ax scenario with following choice 
of the soft masses for the third generation squarks: Mq 3 = M^ 3 = = 800 GeV, and 
X t — A t — n cot (3 = —1600 GeV. The advantage of such a choice is that other SUSY param- 
eters more easily survive the constraints (so that the bounds we obtain are conservative). 
Moreover, we assume the grand unification relation for the gaugino masses: 

Mi : M 2 : M 3 ~ 1 : 1.83 : 5.26 . (17) 

This relation is often assumed in studies of SUSY at the TeV scale for it can be easily 
generated in the mSUGRA model [24]. Note that relaxing this relation will give a large 
effect on the light neutralino scenario [25]. 

The surviving points for the three model are displayed in Fig. 1 for the spin-independent 
elastic cross section of x-nucleon scattering. We see that for each model the CDMS II limits 
can exclude a large part of the parameter space allowed by current collider constraints and 
the future SuperCDMS (25 kg) limits can cover most of the allowed parameter space. For 
the MSSM and the NMSSM dark matter mass is roughly in range of 50-400 GeV, while 
for the nMSSM dark matter mass is constrained below 40 GeV by current experiments and 
further constrained below 20GeV by SuperCDMS in case of non-observation. 

From Fig. 1, we can see that the x-nucleon scattering plot of the MSSM and the NMSSM 
are very similar to each other, but very different from nMSSM. This implies that under 
the experiment constraints, the singlet sector will decouple from the MSSM-like sector in 
the NMSSM, then the NMSSM will perform almost the same as the MSSM, However, the 
singlet components change EW scale phenomenology greatly in the nMSSM. This can also 
be seen in Fig. 2 and Fig. 3. We can see that for both the MSSM and the NMSSM x°i is 
bino-dominant, while for the nMSSM Xi is singlino-dominant, and the region allowed by 
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FIG. 1: The scatter plots (taken for Ref. [12]) for the spin-independent elastic cross section of 
X-nucleon scattering. The '+' points (red) are excluded by CDMS limits (solid line), the 'x' (blue) 
would be further excluded by SuperCDMS 25kg [20] in case of non-observation (dash-dotted line), 
and the 'o' (green) are beyond the SuperCDMS sensitivity. 
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FIG. 2: Same as Fig. 1, but projected on the plane of |-/Vn| 2 and I-/V15I 2 versus dark matter mass, 
(taken for Ref. [12]) 

CDMS limits (and SuperCDMS limits in case of non-observation) favors a more bino-like Xi 
for the MSSM/NMSSM and a more singlmo-like xl for the nMSSM. For the MSSM/NMSSM 
the LSP lower bound around 50 GeV is from the chargino lower bound of 103.5 GeV plus 
the assumed GUT relation Mi ~ O.5M2; while the upper bound around 400 GeV is from the 
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FIG. 3: Same as Fig. 1, but showing the chargino mass m + versus the LSP mass. The dashed 
lines indicate m„ + = m v o. (taken for Ref. [121) 



bino nature of the LSP (Mi cannot be too large, must be much smaller than other relevant 
parameters) plus the experimental constraints like the muon g-2 and B-physics. If we do 
not assume the GUT relation Mi ~ 0.5M2, then Mi can be as small as 40 GeV and the LSP 
lower bound in the MSSM/NMSSM will not be sharply at 50 GeV. (We talk about it in the 
following section.) For both the MSSM and the NMSSM, the CDMS limits tend to favor 
a heavier chargino and ultimately the SuperCDMS limits tend to favor a wino-dominant 
chargino with mass about 2m x o. Note that, there still can be a singlino dominant LSP in 
some parameter space of the NMSSM [26], but in the scan range Eq. (16) listed above, 
getting such singlino dominant LSP needs some fine-tuning, thus we do not focus on it. 

In Fig. 4 we show the value of |A| versus the charged Higgs mass in the NMSSM and 
the nMSSM. This figure indicates that A larger than 0.4 is disfavored by the NMSSM. The 
underlying reason is that hiXiXi depends on A explicitly and large A can enhance x-nucleon 
scattering rate. By contrast, although CDMS has excluded some points with large A in the 
nMSSM, there are still many surviving points with A as large as 0.7. We have talked the 
reason above: to get a heavy Xi, one need a large A, a small tan f3 as well as a moderate n e g. 

From the survived parameter space for all the model above, we should know that the 
Higgs decay will be similar for the MSSM and the NMSSM, but quite different from the 
nMSSM. This can be seen in Fig. 5 which shows decay branching ratio of hi — > XiXi versus 
the mass of the SM-like Higgs boson hsM ( which is hi here, and it is Higgs doublet H u and 
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FIG. 4: Same as Fig. 1, but projected on the plane of |A| versus the charged Higgs mass in the 
NMSSM and the nMSSM. (taken for Ref. [12]) 




m h (GeV) 



m h (GeV) 



m h (GeV) 



FIG. 5: Same as Fig. 1, but projected for the decay branching ratio of hsM — > XiXi versus the 
mass of the Higgs boson /ism- (taken for Ref. [12]) 

Hd dominant ). Such a decay is strongly correlated to the x-nucleon scattering because the 
coupling hiXiXi is involved in both processes. We see that in the MSSM and the NMSSM 
this decay mode can open only in a very narrow parameter space since Xi cannot be so 
light, and in the allowed region this decay has a very small branching ratio (below 10%). 
However, in the nMSSM this decay can open in a large part of the parameter space since 
the LSP can be very light, and its branching ratio can be quite large (over 80% or 90%). 
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B. light dark matter in the NMSSM 



As talked in the introduction, the data of CoGeNT experiment favors a light dark matter 
around 10 GeV. However, we scan the parameter space in the MSSM and find that it is 
very difficult to find a neutralino Xi lighter than about 28 GeV, unless when it is associated 
with a light stau as the next to the lightest supersymmetric particle (NLSP), but such 
scenario always needs a fine-tuning in the parameter space [27]. The main reason for the 
absence of a lighter x° is that the dominant annihilation channel for Xi in the early universe 
is XiXi ~ >* bb through s-channel exchange of the pseudoscalar Higgs boson (A) and the 
measured dark matter relic density requires uia ~ (90-100) GeV and tan/3 ~ 50, this is in 
conflict with the constraints from the LEP experiment and B physics [28-30]. The LHC 
data gives an even more stronger constraint on the light pseudoscalar scenario [31] such that 
light dark matter seems impossible in the MSSM. Though in the nMSSM the neutralino Xi 
can be as light as 10 GeV (shown in Fig. 1), the scattering rate is much lower under the 
CoGeNT-favored region. In the NMSSM, however, with the participation of singlet sector 
one can get very light [4] Higgs. This feature is particularly useful for light Xi scenario since 
it opens up new important annihilation channels for xi, i-e., either into a pair of hi (or a±) or 
into a pair of fermions via s-channel exchange of hi (or ai) [30, 32, 33]. For the former case, 
Xi must be heavier than hi (ai); while for the latter case, due to the very weak couplings of 
hi (ai) with Xi an d with the SM fermions, a resonance enhancement (i.e. m/ ll or m ai must 
be close to 2m^o) is needed to accelerate the annihilation. So a light Xi ma y be necessarily 
accompanied by a light hi or ai to provide the required dark matter relic density. From the 
discussion in the upper section, light Xi can be obtained by releasing the GUT relation Eq. 
(17), thus LSP in the NMSSM may explain the detection of CoGeNT. Note that, as the 
LSP in the nMSSM is singlino dominant, relaxing the GUT relation will not the change the 
phenomenology of dark matter and Higgs too much. 

Now we discuss how to get a light hi or a± in the NMSSM. A light a\ can be easily 
obtained when the theory is close to the U(l)# or U(1)pq symmetry limit, which can be 
realized by setting the product kA k to be negatively small [4]. In contrast, a light hi can not 
be obtained easily. However, as shown below, it can still be achieved by somewhat subtle 
cancelation via tuning the value of A K . We note that for any theory with multiple Higgs 
fields, the existence of a massless Higgs boson implies the vanishing of the determinant of 
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FIG. 6: The scatter plots (taken for Ref. [13]) of the parameter samples which survive all con- 
straints, with 'x' (red) and 'A' (green) corresponding to a light h\ and a light a±, respectively. 
The left frame is projected on the <r SI -m x plane, while the right frame is projected on the a-m^ 
plane (denoted by 'x') and the a -m ai plane (denoted by 'A'). The curves are the limits from 
CoGeNT [9], CDMS [7], while the contour is the CoGeNT-favored region [9]. 

its squared mass matrix and vice versa. For the NMSSM, at tree level the parameter A K 
only enters the mass term of the singlet Higgs bosons, so the determinant (Det M 2 ) of the 
mass matrix of the CP-even Higgs bosons depends on A K linearly [4]. When other relevant 
parameters are fixed, one can then obtain a light hi by varying A K around the value A K which 
is the solution to the equation Det.M 2 = 0. In practice, one must include the important 
radiative corrections to the Higgs mass matrix, which will complicate the dependence of Ai 2 
on A K . However, we checked that the linear dependence is approximately maintained by 
choosing the other relevant parameters at the SUSY scale, and one can solve the equation 
iteratively to get the solution A K . 

In Fig. 6 we display the surviving parameter samples, showing the x-nucleon scattering 
cross section versus the neutralino dark matter mass (left frame) and versus the mass of 
hi or ai (right frame). It shows that the scattering rate of the light dark matter can reach 
the sensitivity of CDMS and, consequently, a sizable parameter space is excluded by the 
CDMS data [20] . The future CDMS experiment can further explore (but cannot completely 
cover) the remained parameter space. Note that in the light- hi case the scattering rate 
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can be large enough to reach the sensitivity of CoGeNT and to cover the CoGeNT-favored 
region. The underlying reason is that the x-nucleon scattering can proceed through the 
t-channel exchange of the CP-even Higgs bosons, which can be enhanced by a factor 
for a light hi [32]; while a light ai can not give such an enhancement because the CP-odd 
Higgs bosons do not contribute to the scattering in this way. We noticed that the studies in 
[29, 34] claimed that the NMSSM is unable to explain the CoGeNT data because they did 
not consider the light-/^ case. 

In the light x\ scenario, h$M m &y decay exotically into h\hi or a±ai, and conse- 

quently the conventional decays are reduced. This feature is illustrated in Fig. 7, which 
shows that the sum of the exotic decay branching ratios may exceed 50% and the traditional 
decays h$ M — > 66, tt, WW*, 77 can be severely suppressed. Numerically, we find that the 
branching ratio of h$M — > bb is suppressed to be below 30% for all the surviving samples in 
the light- /ii (h 2 is h SM ) case and for about 96% of the surviving samples in the light-ai (hi 
is ^sm) case (for the remaining 4% of the surviving samples in the light-ai case, the decay 
hsM — > oiOi is usually kinematically forbidden so that the ratio of hsyi — > bb may exceed 
60%). Another interesting feature shown in Fig. 7 is that, due to the open-up of the exotic 
decays, h$M ma y be significantly lighter than the LEP bound. This situation is favored by 
the fit of the precision electro- weak data and is of great theoretical interest [35] . 

Since the conventional decay modes of h SM may be greatly suppressed, especially in the 
light-/ii case which can give a rather large x-nucleon scattering rate, the LHC search for hsM 
via the traditional channels may become difficult. Now the LHC observed a new particle 
in the mass region around 125-126 GeV which is the most probable the long sought Higgs 
boson [36]. In this mass range, the most important discovering channel of h SM at the LHC 
is the di-photon signal. In Fig. 8 we give the ratio of the di-photon production rate to 
the SM at the LHC with yfs = 7 TeV. In calculating the rate, we used the narrow width 
approximation and only considered the leading contributions to pp — > h$M from top quark, 
bottom quark and the squark loops. 

Fig. 8 indicates that, compared with the SM prediction, the ratio in the NMSSM in the 
light Xi scenario is suppressed to be less than 0.4 for the light-/ii case. For the light-ai case, 
most samples (about 96%) predict the same conclusion. Since in the light-/ii case the x~ 
nucleon scattering rate can reach the CoGeNT sensitivity, this means that in the framework 
of the NMSSM the CoGeNT search for the light dark matter will be correlated with the 
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FIG. 7: Same as Fig. 6, but showing the decay branching ratios of the SM-like Higgs boson h^M- 

Here Br(hsM — > XiX°j) denotes the total rates for all possible hsM xfx? decays, (taken for Ref. 



[13]) 



LHC search for the Higgs boson via the di-photon channel. We checked that, once the future 
XENON experiment fails in observing dark matter, less than 1% of the surviving samples 
in light <2i case predict the ratio of di-photon signal larger than 0.4. 
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h SM mass 

FIG. 8: Same as Fig. 6, but showing the diphoton production rate of the SM-like Higgs boson at 
the LHC. 

IV. GENERAL EXTENSION FOR THE EXPLANATION TO PAMELA 

To explain the PAMELA excess by dark matter annihilation, there are some challenges. 
First, dark matter must annihilate dominantly into leptons since PAMELA has observed no 
excess of anti-protons [11] (However, as pointed in Ref. [37], this statement may be not so 
solid due to the significant astrophysical uncertainties associated with their propagation). 
Second, the explanation of PAMELA excess requires an annihilation rate which is too large 
to explain the relic abundance if dark matter is produced thermally in the early universe. 
To tackle these difficulties, a new theory of dark matter was proposed in Ref. [38]. In this 
new theory the Sommerfeld effect of a new force in the dark sector can greatly enhance 
the annihilation rate when the velocity of dark matter is much smaller than the velocity at 
freeze-out in the early universe, and dark matter annihilates into light particles which are 
kinematically allowed to decay to muons or electrons. 

The above fancy idea is hard to realize in the MSSM, because there is not a new force 
in the neutralino dark matter sector to induce the Sommerfeld enhancement and neutralino 
dark matter annihilates largely to final states consisting of heavy quarks or gauge and/or 
Higgs bosons [19, 39]. However, as discussed in Ref. [40], in a general extension of the 
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MSSM by introducing a singlet Higgs superfield, the idea in Ref. [38] can be realized by the 
singlino-like neutralino dark matter: 

(i) The singlino dark matter annihilates to the light singlet Higgs bosons and the relic 
density can be naturally obtained from the interaction between singlino and singlet 
Higgs bosons. 

(ii) The singlet Higgs bosons, not related to electro-weak symmetry breaking, can be 
light enough to be kinematically allowed to decay dominantly into muons or electrons 
through the tiny mixing with the Higgs doublets. 

(iii) The Sommerfeld enhancement needed in dark matter annihilation for the explanation 
of PAMELA result can be induced by the light singlet Higgs boson. 

In the following section, we will show how does this happen, the Higgs decay are also 
investigated. 

A. Higgs and neutralinos spectrum 

If introduce a singlet Higgs to the MSSM in general, the renormalizable holomorphic 
superpotential of Higgs is given by Ref. [40] 

W = f iH u -H d + \SH u -H d + r ] S+^i2 s S 2 + ^KS 3 , (18) 

which include linear term, quadratic term, cubic term of singlet superfield (like Wess-Zumino 
model [41]). Note that in such case, we do not require the singlet to solve the \x problem. 
The soft SUSY-breaking terms are given by 

V soit = m 2 u \H u \ 2 + m 2 d \H d \ 2 + m 2 s \S\ 2 

+ (BfiH u ■ H d + \A\ H u ■ H d S + C V S+\ B sf i s S 2 + \ kA k S 3 + h.c.) . (19) 

After the Higgs fields develop the VEVs v u , v d and s, i.e., we get the similar Higgs spectrum 
as the NMSSM and the nMSSM which is 
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(1) The CP-even Higgs mass matrix in the basis (</>„, <p d) c) is given by 

M h>u = g 2 v 2 u + cot (3 [Xs(A x + ks + // s ) + B/x] , (20) 
M K 22 = g 2 v\ + tan p [Xs(A x + ks + /i s ) + £//] , (21) 
M h ,ss = X(A X + ^)— " A-(«S + «2) + + 4^ + 3/i s ) - ^, (22) 

5 5 5 

= (2A 2 - g 2 )v u v d - Xs{A x + ks + n,) - Bp, (23) 
M h ,i3 = 2X(fx + Xs)v u - Xv d (A x + 2ns + fx s ), (24) 
M h ,23 = 2X(n + Xs)v d - Xv u (A x + 2ks + fi s ), (25) 

where g 2 = (gf + g 2 )/2 with g 1 and g 2 being respectively the coupling constant of 
SU(2) and U(l) in the SM. 

(2) The CP-odd Higgs mass matrix M a is given by 

M a ,u = (tan (3 + cot (3)[Xs(A x + ks + fi s ) + Bfj], (26) 

M a ,22 = 4XKv u v d + X(A x + /2 s )^-X% 2 u + v 2 d ) 

Cn 

-ks(3A k + fi s ) '-- 2B s fi s , (27) 

s 

M a ,12 = X^vl+v 2 d (A x -2Ks-ii s ). (28) 

Note that here we have dropped the Goldstone mode, thus there left a 2 x 2 mass 
matrix in the basis (A,£). and it can be diagonalized by an orthogonal 2x2 matrix 
P' and the physical CP-odd states are given by (ordered as m ai < m a2 ) 

a x = P' U A + P[ 2 S: = P^cos ^ u + sin p<p d ) + P[£, (29) 
a 2 = P' 21 A + P' 22 Sj = P^(cos/3^ + sin/3^ d ) + P^. (30) 

(3) The charged Higgs mass matrix M± in the basis H d ) is given by 

M±= (xs(A x + Ks + fi s ) + Bfi + h u h d $ -X 2 )) l COt(3 1 |, (31) 
V 2 / \ 1 tan/3 
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(4) The neutralino mass matrix is : 

/ Mi 



Mr 



m z s w s b -m z s w cb 

M 2 -m z c w s h m z c w c b 

m z swSb —iTL Z swSb —fj, 

—mzSwCb —m z cwCb —[i 






—Xvcb 
—Xvsb 



\ 











-Xvcb 



(32) 



-Xvsb 2ks + ji s J 



B. Explanation of PAMELA and implication on Higgs decays 

To explain the observation of PAMELA, a± is singlet-dominant, while hi is singlet- 
dominant and the next-to-lightest h 2 is doublet-dominant (/ism)- We use the notation: 



a = ai, h = hi, h. 



SM 



ho. 



(33) 



As discussed in Ref. [40], when the lightest neutralino Xi m Eq. (11) is singlino-dominant, 
it can be a perfect candidate for dark matter. As shown in Fig. 9, such singlino dark 
matter annihilates to a pair of light singlet Higgs bosons followed by the decay h — >■ aa (h 
has very small mixing with the Higgs doublets and thus has very small couplings to the 
SM fermions). In order to decay dominantly into muons, a must be light enough. Further, 
in order to induce the Sommerfeld enhancement, h must also be light enough. From the 
superpotential term kS 3 we know that the couplings hxiXi an d <^XiXi are proportional to k. 
To obtain the relic density of dark matter, k should be 0(1). h, a are singlet- dominant and 
Xi is singlino-dominant, this implies small mixing between singlet and doublet Higgs fields. 
From the superpotential in Eq.(18) we see that this means the mixing parameter A must be 
small enough. On the other hand, from Eq. (22) and Eq. (27) lightness of hi and ai also 
require A and other term approaching to zero. Therefore, in our scan we require parameters 
A K and B s has the relation: 



A K ~ ( -Aks - 3/jl s + ^\ ) , 

KS 



2B s [i s 



1A Cl l I 

s 



(34) 
(35) 



to realize light hi and a x . 

The numerical results of this model are displayed in different planes in Figs. 10-12. We 
see from Fig. 10 that in the range 2m ^ < m a < 2m n , a decays dominantly into muons. It 
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FIG. 9: Feynman diagrams for singlino dark matter annihilation where Sommerfeld enhancement 
is induced by exchanging h. (taken from Ref. [14]) 
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m a (GeV) 

FIG. 10: The scatter plots showing the decay branching ratios a — > n + fJ-~ (muon), a — >• gg (gluon) 
and a — >■ ss (s-quark) versus m a for A = 10 -3 . (taken from Ref. [14]) 

is clear that h can be as light as a few GeV, which is light enough to induce the necessary 
Sommerfeld enhancement as shown in Fig. 11. 

In left plot of Fig. 12, we show the branching ratios of hsM decays. We see that in 
the allowed parameter space hsu tends to decay into aa or hh instead of bb. This can 
be understood as following, the MSSM parameter space is stringently constrained by the 
LEP experiments if hsu is relatively light and decays dominantly to bb, and to escape such 
stringent constraints hsM tends to have exotic decays into aa or hh. As a result, the allowed 
parameter space tends to favor a large A K , as shown in right plot of Fig. 12, which greatly 
enhances the couplings hsuo-a and h^uhh through the soft term kA k S 3 although S has a 
small mixing with the doublet Higgs bosons. Such an enhancement can be easily seen. Take 
the coupling ho,y[hh as an example, the soft term kA k S 3 gives a term nA K a 3 which then 
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m h (GeV) 

FIG. 11: Same as Fig. 10, but showing the Sommerfeld enhancement factor induced by h. (taken 
from Ref. [14]) 




FIG. 12: Same as Fig. 10, but showing branching ratio of /ism aa,hh versus mh 3m and \A K \ 
versus the branching ratio of hsM — > aa, hh. (taken from Ref. [14]) 

gives the interaction kA k {U^) 2 U^ h^hh because a = Uf^hi + U^h 2 + U^h 3 with h\ = h 
and h>2 = hgM (see Eqs. (9) and (33)). Although the mixing (U^) 2 U^ is small for a small 
A, a large A K can enhance the coupling h^uhh. Note that as the mass of the observed Higgs 
boson at the LHC is around 125 GeV, thus in the MSSM, the dominant decay mode of 
h§M is bb. In this general singlet extension of the MSSM, its dominant decay mode may be 
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changed to aa or hh, as shown in our above results. 

Finally, we note that for the specified singlet extensions like the nMSSM and the NMSSM, 
the explanation of PAMELA and relic density through Sommerfeld enhancement is not 
possible. The reason is that the parameter space of such models is stringently constrained 
by various experiments and dark matter relic density as shown in the above section, and, 
as a result, the neutralino dark matter may explain either the relic density or PAMELA, 
but impossible to explain both via Sommerfeld enhancement. For example, in the nMSSM 
various experiments and dark matter relic density constrain the neutralino dark matter 
particle in a narrow mass range [42], which is too light to explain PAMELA. 

V. SUMMARY 

At last we summarize here, the SUSY dark matter and Higgs physics will be changed if 
introducing a singlet to the MSSM. Under the latest results of dark matter detection, we 
have: 

1. In the MSSM, the NMSSM and the nMSSM, the latest detection result can exclude 
a large part of the parameter space allowed by current collider constraints and the 
future SuperCDMS and XENON can cover most of the allowed parameter space. 

2. Under the new dark matter constraints, the singlet sector will decouple from the 
MSSM-like sector in the NMSSM, thus the phenomenologies of dark matter and Higgs 
are similar to the MSSM. The singlet sector make the nMSSM quite different from 
the MSSM, the LSP in the nMSSM are singlet dominant, and the SM-like Higgs will 
mainly decay into the singlet sector. Future precision measurements will give us an 
opportunity to determine whether the new scalar is from standard model or from 
SUSY. Perhaps the nMSSM will be the first model be excluded for its much larger 
branching ratio of invisible Higgs decay. 

3. The NMSSM can allow light dark matter at several GeV exists. Light CP-even or 
CP-odd Higgs boson must be present so as to satisfy the measured dark matter relic 
density. In case of the presence of a light CP-even Higgs boson, the light neutralino 
dark matter can explain the CoGeNT and DAMA/LIBRA results. Further, we find 
that in such a scenario the SM-like Higgs boson will decay predominantly into a pair 
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of light Higgs bosons or a pair of neutralinos and the conventional decay modes will 
be greatly suppressed. 

4. The general singlet extension of the MSSM gives a perfect explanation for both the 
relic density and the PAMELA result through the Sommerfeld enhanced annihilation 
into singlet Higgs bosons (a or h followed by h — > aa) with a being light enough to 
decay dominantly to muons or electrons. Although the light singlet Higgs bosons have 
small mixing with the Higgs doublets in the allowed parameter space, their couplings 
with the SM-like Higgs boson hsu can be enhanced by the soft parameter A K . In 
order to meet the stringent LEP constraints, the Hsm tends to decay into the singlet 
Higgs pairs aa or hh instead of bb. 
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